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We present the first observation of the Z'y — > vv^ process at the Tevatron at 5.1 standard 
deviations significance, based on 3.6 fb -1 of integrated luminosity collected with the DO detector at 
the Fermilab Tevatron pp Collider at y^s — 1.96 TeV. The measured Zy cross section multiplied by 
the branching fraction of Z — » vv is 32 ± 9(stat. + syst.) ± 2(lumi.) fb for the photon Et > 90 GeV. 
It is in agreement with the standard model prediction of 39 ±4 fb. We set the most restrictive limits 
on anomalous trilinear Z77 and ZZj gauge boson couplings at a hadron collider to date, with three 
constraints being the world's strongest. 



PACS numbers: 12.15.Ji, 12.60.Cn, 13.38.Dg, 13.40.Em, 13.85.Qk, 14.70.-e 



The standard model (SM) of electrowcak interactions U(l). The symmetry transformations of the group allow 
is described by the non-Abelian gauge group SU (2) x interactions involving three gauge bosons (7, W, and Z) 
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through trilinear gauge boson couplings. However, the 
SM forbids such vertices for the photon and the Z boson 
at the lowest "tree" level, i.e., the values of the Z77 and 
ZZ7 couplings vanish. The cross section for the SM Z7 
production is very small. However, the presence of finite 
(anomalous) Zjj and ZZj couplings can enhance the 
yields, especially at higher values of the photon trans- 
verse energy (Et). As we are marginally sensitive to 
one-loop SM contributions [3, 0] to these vertices, obser- 
vation of an anomalously high Z7 production rate could, 
therefore, indicate the presence of new physics. 

To preserve S'-matrix unitarity, the anomalous cou- 
plings must vanish at high center-of-mass energies. 
Hence, the dependence on the center-of-mass energy 
has to be included in the definition of such couplings. 
This can be done by using a set of eight complex pa- 
rameters hj(i = 1, ...,4;V = Z, 7) of the form hj = 
hY Q /(l + s/A 2 ) n Here, s is the square of the center-of- 
mass energy in the partonic collision, A is a scale related 
to the mass of the new physics responsible for anomalous 
Zj production, and hY is the low energy approximation 
of the coupling. Following Ref. we will use n = 3 for 
h\ and h\ , and n — 4 for hX and h\ ■ This choice of n 
guarantees the preservation of partial-wave unitarity, and 
makes the vertex function terms proportional to h\ and 

behave in the same way as terms proportional to 
and h\ at high energies. Couplings h\ and h^o (h^o and 
hYo) are CP-violating (CP-conserving). In this Letter, 
we set limits on the size of the real parts of the anoma- 
lous couplings: Re(hY ), which we refer to as ATGC in 
the following. 

In the past, studies of Z7 production have been per- 
formed by the CDF 4] and DO 0, 0] collaborations at 
the Tevatron Collider, as well as at the CERN LEP Col- 
lider by the L3 7], and OPAL (§] collaborations. The 
most recent combination of LEP results can be found in 
Ref. 0. 

The DO detector [l(| consists of a central-tracking sys- 
tem, liquid-argon/uranium calorimeters, and a muon sys- 
tem. The tracking system comprises a silicon microstrip 
tracker (SMT) and a central fiber tracker (CFT), both lo- 
cated within a « 2 T superconducting solenoid, and pro- 
vides tracking and vertexing up to pseudorapidities [11 1 
of \n\ w 3.0 and |?7| « 2.5, respectively. The central 
and forward preshower detectors (CPS and FPS) are lo- 
cated between the superconducting coil and the calorime- 
ters, and consist of three and four layers of scintillator 
strips, respectively. The liquid- argon/uranium calorime- 
ter is divided into a central calorimeter (CC) and two 
end calorimeters (EC), covering pseudorapidities up to 
\n\ Kt 1.1 and \r)\ w 4.2, respectively. The calorimeters are 
segmented into an electromagnetic section (EM), com- 
prised of four layers, and a hadronic section, divided lon- 
gitudinally into fine and coarse sections. The calorimeter 
is followed by the muon system, consisting of three lay- 
ers of tracking detectors and scintillation trigger counters 



and a 1.8 T iron toroidal magnet located between the two 
innermost layers. The muon system provides coverage to 
1 77 1 2. Arrays of plastic scintillators in front of the EC 
cryostats are used to measure the luminosity. 

Data for this analysis were collected with the DO de- 
tector in the period from 2002 to 2008, and correspond 
to an integrated luminosity of 3.6 fb _1 after the appli- 
cation of data-quality and trigger requirements. Events 
must satisfy a trigger from a set of high-E'r single EM- 
cluster triggers, which are (99 ± 1)% efficient for photons 
of E T > 90 GeV. 

Photons are identified as calorimeter clusters with at 
least 95% of their energy deposited in the EM calorime- 
ter, with transverse and longitudinal distributions con- 
sistent with those of a photon, and spatially isolated 
in the calorimeter and in the tracker. A cluster is iso- 
lated in the calorimeter if the isolation variable 2 = 
[£ tot (0.4)-£ E M(0.2)]/£ EM (0.2) < 0.07. Here, E tot (0.4) 
is the total energy (corrected for the contribution from 
multiple pp interactions) deposited in a calorimeter cone 
of radius TZ = ^/(Ary) 2 + (A0) 2 = 0.4, and £ E m(0.2) is 
the EM energy in a cone of radius TZ — 0.2. The track 
isolation variable, defined as the scalar sum of the trans- 
verse momenta of all tracks that originate from the in- 
teraction vertex in an annulus of 0.05 < TZ < 0.4 around 
the cluster, must be less than 2 GeV. 

We obtain the photon sample by selecting events with 
a single photon candidate of Et > 90 GeV and |r/| < 1.1, 
and require a missing transverse energy in the event of 
$ T > 70 GeV, which effectively suppresses the multijet 
background. The Ifi T is computed as the negative vector 
sum of the Et of calorimeter cells and corrected for the 
transverse momentum of reconstructed muons and the 
energy corrections to reconstructed electrons and jets. 
To minimize large Tft T from mismeasurement of jet en- 
ergy, we reject events with jets with Et > 15 GeV. 
We also reject events containing reconstructed muons, 
and events with cosmic-ray muons identified through a 
timing of their signal in the muon scintillators. Events 
with additional EM objects with Et > 15 GeV are 
rejected. To suppress W boson decays into leptons, 
events with reconstructed high-p^ tracks are removed. 
To reduce the copious non-collision background (events 
in which muons from the beam halo or cosmic rays un- 
dergo bremsstrahlung, and produce energetic photons), 
we use a pointing algorithm 12J, exploiting the trans- 
verse and longitudinal energy distribution in the EM 
calorimeter and CPS. This algorithm is based on esti- 
mates of z positions of production vertices (zem) along 
the beam direction assuming that given EM showers are 
initiated by photons, and utilizes the distance of closest 
approach (DCA) [13j of the direction of the EM shower 
to the z axis. We require |zem — zsj\ < 10 cm, where 
zem is the z position of the interaction vertex predicted 
by the pointing algorithm and zy is the z position of the 
chosen (often nearest) reconstructed vertex. 
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TABLE I: Summary of background estimates, and the number 
of observed and SM predicted events. 



Number of events 


W —* ev 


9.67 ± 0.30(stat.) ± 0.48(syst.) 


non- collision 


5.33 ± 0.39(stat.) ± 1.91(syst.) 


W/Z + jet 


1.37 ± 0.26(stat.) ± 0.91(syst.) 




0.90 ± 0.07(stat.) ± 0.12(syst.) 


Total background 


17.3 ± 0.6(stat.) ± 2.3(syst.) 


N SM 


33.7 ±3.4 


N obs 


51 



Following the procedure described in Ref. [14j |. we 
estimate the fraction of non-collision and W/Z events 
with misidcntified jets backgrounds in the final candidate 
events by fitting their DCA distribution to a linear sum 
of three DCA templates. These templates are: a tem- 
plate resembling the signal, a non-collision template, and 
a misidentified jets template. Most of the signal photons 
are concentrated in the region with DCA < 4 cm. There- 
fore, we restrict the analysis to this particular range. 

Other backgrounds to the 7+^ T signal arise from clec- 
troweak processes such as W — ► ev, where the electron is 
misidentified as a photon due to inefficiency of the tracker 
or hard bremsstrahlung, and W"f, where the lepton from 
the W boson decay is not reconstructed. 

The W — > ev background is estimated using a sam- 
ple of isolated electrons. We apply the same kinematic 
requirements as in the photon sample, and scale the 
remaining number of events by the measured rate of 
electron-photon misidentification, which is 0.014 ±0.001. 
The Wj background is estimated using a sample of 
Monte Carlo (MC) events generated with pythia [la| . 
These events are passed through a detector simulation 
chain based on the GEANT package and recon- 

structed using the same software as used for data. Af- 
ter imposing the same selection requirements as for the 
photon sample, scale factors are applied to correct for 
differences between simulation and data. The summary 
of backgrounds is shown in Table fl] 

After applying all selection criteria, we observe 
51 candidate events with a predicted background of 
17.3 ± 0.6(stat.) ± 2.3(syst.) events. To estimate the to- 
tal acceptance of the event selection requirements, we use 
MC samples produced with a leading-order (LO) gen- 
erator [3|] , passed through a parameterized simulation of 
the DO detector. The next-to-leading order (NLO) QCD 
corrections arising from soft gluon radiation and virtual 
one-loop corrections are taken into account through the 
adjustment of the photon Et spectrum using a AT-factor, 
estimated using a NLO Z7 event generator 17J. As we 
require no jets with Et > 15 GeV to be present in the 
final state, the NLO corrections, integrated over the pho- 
ton Et range after the photon Et > 90 GeV requirement, 
are « 2% or smaller both for the SM and anomalous Zj 



production. The NLO corrections distribution is fitted 
with a smooth function, with an uncertainty of ~ 5% 
arising from the fit. The uncertainty on the if- factor 
from the jet energy scale and resolution is estimated to 
be « 3%. Based on this simulation, the expected number 
of events from the SM signal is estimated to be 33.7 ± 3.4 
events. The number of observed events (A^bs) and the 
number of predicted events (N„^~.) are summarized in 
Table H 

The Z"f cross section multiplied by the branching frac- 
tion of Z — > vv is measured to be 32 ± 9(stat. ± syst.) ± 
2(lumi.) fb for the photon E T > 90 GeV, which 
is in goo d ag reement with the NLO cross section of 
39 ± 4 fb 17]. The main contribution to the total un- 
certainty on the measured cross section is the statistical 
uncertainty on the small number of events in the final 
sample, and is a factor of four to five larger than the 
individual systematic uncertainties on photon identifi- 
cation, choice of parton distribution functions (PDF), 
and kinematic criteria. The uncertainty on the theo- 
retical cross section comes mainly from the choice of 
PDF (7%) and estimation of the NLO A'-factor (5.5%). 
To estimate the statistical significance of the measured 
cross section, we perform 10 8 background-only pseudo- 
experiments and calculate the p-value as the fraction 
of pseudo-experiments with an estimated cross section 
above the measured one. This probability is found to be 
3.1 x 10 -7 , which corresponds to a statistical significance 
of 5.1 standard deviations (s.d.), making this the first 
observation of the Z7 — > vv^j process at the Tevatron. 

To set limits on the ATGC, we compare the photon 
Et spectrum in data with that from the sum of expected 
Z^f signal 0,[l3| and the background (see Fig. [I]) for each 
pair of couplings for a grid in which h^ runs from -0.12 
to 0.12 with a step of 0.01, and h\ varies from -0.08 to 
0.08 with a step of 0.001. The MC samples are generated 
with the LO Z7 generator (corrected for the NLO effects 
with an independent A'-factor for the form-factor 
A = 1.5TeV. 

Assuming Poisson statistics for the signal and Gaus- 
sian distribution of all the systematic uncertainties on 
the generated samples and on the backgrounds, we cal- 
culate the likelihood of the photon Et distribution in 
data given the prediction for hypothesized ATGC. To 
set limits on any individual ATGC at the 95% confi- 
dence level (C.L.), we set the other anomalous couplings 
to zero. The resulting limits in the neutrino channel 
alone are \hj \ < 0.036, \hj \ < 0.0019 and |/if | < 0.035, 
\hf \ < 0.0019. To further improve the sensitivity, we 
generate the Z-f — » ££-f (£ = e, fi) MC samples for 
these couplings and A = 1.5 TeV, and set limits on 
ATGC for the 1 fb -1 data sample used in the previous 
Z7 analysis The combination of all three channels 
yields the most stringent limits on the ATGC set at a 
hadron collider to date: \hj \ < 0.033, \KU < 0.0017 and 



|/if | < 0.033, 



K Z ol 



< 0.0017. This is roughly a factor of 
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three improvement over the results published in Ref . [6j . 
The limits on the hf , hf , and hj couplings improve 
on the constraints from LEP2, and are the most restric- 
tive to date. The limits on the CP-violating couplings h\ 
and hYa are, within the precision of this measurement, the 
same as the limits on h^ and h^ , respectively. Hence, 
we can constrain the strength of the couplings but not 
the phase. As the described method is sensitive only to 
the magnitude and the relative sign between couplings, 
the one- and two-dimensional limits are symmetric with 
respect to the SM coupling under simultaneous exchange 
of all signs. The 95% C.L. one-dimensional limits and 
two-dimensional contours are shown in Figs. [5Ji and [2b 
for the CP-conserving Z77 and ZZy couplings, respec- 
tively. 
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FIG. 1: Photon Et spectrum in data (solid circles), sum of 
backgrounds (dash-dot line), and sum of MC signal and back- 
ground for the SM prediction (solid line) and for the ATGC 
prediction with hj = 0.09 and hj = 0.005 (dashed line). 
The shaded band corresponds to the ± 1 s.d. total uncer- 
tainty on the predicted sum of SM signal and background. 



In summary, we observe 51 vvj candidates with 
17.3 ± 0.6(stat.) ± 2.3(syst.) background events us- 
ing 3.6 fb -1 of data collected with the DO detector at 
the Tevatron. We measure the most precise Z7 — > 
vv^j cross section to date at a hadron collider of 32 ± 
9(stat. + syst.)±2(lumi.) fb for the photon E T > 90 GeV, 



in agreement with the SM prediction of 39 ± 4 fb [17 1. 
The statistical significance of this measurement is 5.1 s.d., 
making it the first observation of the Z"f — > vvy process 
at the Tevatron. We set the most restrictive limits on 
the real parts of the anomalous trilinear gauge couplings 
at hadron colliders at the 95% C.L. of | 1 < 0.033, 
\hj \ < 0.0017 and \h$ Q \ < 0.033, \hf Q \ < 0.0017. Three 
of these limits are world's best to date. These limits ap- 
proach the range of expectations for the contributions 
due to one-loop diagrams in the SM [l|, Q ■ 
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FIG. 2: Two-dimensional bounds (ellipses) at 95% C.L. on 
CP-conserving (a) Z77 and (b) ZZy couplings. The crosses 
represent the one-dimensional bounds at the 95% C.L. setting 
all other couplings to zero. The dashed lines indicate the 
unitarity limits for A = 1.5 TeV. 
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